Absence of the non-percolating phase for percolation on the non-planar Hanoi network 
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We investigate bond percolation on the non-planar Hanoi network (HN-NP), which was studied in 
[Boettcher et al. Phys. Rev. E 80 (2009) 041115]. We calculate the fractal exponent of a subgraph 
of the HN-NP, which gives a lower bound for the fractal exponent of the original graph. This lower 
bound leads to the conclusion that the original system does not have a non-percolating phase, where 
only finite size clusters exist, for p > 0, or equivalently, that the system exhibits either in the critical 
phase, where infinitely many infinite clusters exist, or the percolating phase, where a unique giant 
component exists. Monte Carlo simulations support our conjecture. 

PACS numbers: 89.75.Hc 64.60.aq 89.65.-s 



I. INTRODUCTION 

Percolation is the simplest model exhibiting a phase 
transition [l[ . Many results for percolation on Euclidean 
lattices have been reported. It is well known that bond 
percolation with open bond probability p on the d{> 2)- 
dimensional Euclidean lattice shows a second order tran- 
sition between the non-percolating phase, where only fi- 
nite size clusters exist, and the percolating phase, where a 
unique giant component almost surely exists, at a unique 
critical point pc- However, this may not be the case for 
non-Euclidean lattices. 

Complex networks have been actively studied in re- 
cent years Among extensive researches carried 
out on complex networks, percolation on various net- 
works has played an important role in clarifying the in- 
terplay between network topology and critical phenom- 
ena [5[. Percolation on uncorrelated networks (repre- 
sented by the configuration model Q) is well described 
by the local tree approximation; there is a phase tran- 
sition between the non-percolating phase and the perco- 
lating phase just as in Euclidean lattice systems, but its 
critical exponents depend crucially on the heterogene- 
ity of the degree distribution of the network 0. On 
the other hand, several authors [8l-IT2| have reported 
that percolation on networks constructed with certain 
growth rules exhibits quite a different phase transition 
from that of uncorrelated networks and Euclidean lat- 
tices, referred to as an infinite order transition with in- 
verted Berezinskii-Kosterlitz-Thouless (BKT) singular- 
ity: (i) The singularity of the phase transition is infinitely 
weak. When p lies above the transition point Pc, the 
order parameter m(p) = limN^oo Snia.xiN;p)/N , where 
Suia.x{N;p) is the mean size of the largest cluster over 



'Electronic address: | hasegawa@m . tohoku. ac.j p | 



percolation trials in the system with N nodes, obeys 
m{p) cx exp[— const. /(Ap)^ ], where Ap = p — Pc- (n) 
Below the transition point, the mean number of clus- 
ters with size s per node obeys the power law, Ug oc s~'^ . 
Furthermore, recent study shows that a hierarchical 
small-world network exhibits a discontinuous transition 
instead of an infinite order transition. 

In [12, [3, [13, it has been found that, for the grow- 
ing network models the region below pc corresponds to 
the critical phase (the intermediate phase), which has 
been observed in nonamenable graphs (NAGs) [H, ■ 
NAGs are defined to be transitive graphs with a posi- 
tive Cheeger constant. Percolation on NAGs (with one 
end) exhibits the following three phases depending on 
the value of p: the non-percolating phase (0 < p < Pci), 
the critical phase (pci < p < Pci), where infinitely 
many infinite clusters exist, and the percolating phase 
(Pc2 < P < !)■ Here, an infinite cluster is defined to be 
a cluster whose size is of order 0{N"-) (0 < a < 1). It 
is called a giant component when a = \. In the critical 
phase, where < a < 1, the system is always in a critical 
state where Ug satisfies a power law [3 ■ 

All of previous studies fsl-fT^. [H, [l^ of percolation on 
growing networks and hierarchical small-world networks 
indicate = Pci < Pc2 < 1 , except in the following case. 
Boettcher et al. investigated bond percolation on the 
non-planar Hanoi network (HN-NP) using the renormal- 
ization group technique [l^. They conclude that there 
are two critical probabilities Pci and Pc2 between zero 
and one: < pd < Pc2 < 1- 

In this paper, we reconsider this model. We show an- 
alytically that the fractal exponent of a subgraph, which 
is a lower bound for that of the HN-NP, takes a non- 
zero value at all p{^ 0), indicating that pd = 0. This 
means that the system is either in the critical phase or 
the percolating phase, not in the non-percolating phase, 
in contrast to the result of [l^. The Monte Carlo simu- 
lations support our analytical prediction. 
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FIG. 1: HN-NP with L = 4 generations. The black-dashed 
lines are the backbones, and the (red-, green-, and blue-)solid 
lines are the shortcuts (of the skeletons Ta(4), Tb{4), and 
Tc(4), respectively). 



II. MODEL 

The HN-NP consists of a one-dimensional chain and 
long-range edges. The HN-NP with L{>2) generation is 
constructed as follows [l^ . (i) Consider a chain of Nl = 
2'^ + 1 nodes. Here, each node 0,1,2, Nl - I) 
connects to node i + We call these edges the backbones. 
(ii) For each combination of i{= 0, 1, 2, • • • ,L ~ 2) and 
j(= 0, 1, 2, • • • , 2-^-^-2 - 1), nodes (4j)2' and (4j + 1)2^ 
are connected to (4j -I- 3)2' and (4j -I- 4)2', respectively. 
We call these edges the shortcuts. The schematic of the 
HN-NP with L = 4 generations is shown in FigH] At 
generation L, the number of backbones is 2^ and the 
number of shortcuts is 2^ — 2 (the total number of edges 
is El is El = 2^+^ — 2). The geometrical properties of 
the HN-NP are as follows [l^ : (i) the degree distribution 
Pk decays exponentially asp2m-i-3 esc 2~™, (ii) the average 
degree (fc) is (fc) = 2El/Nl = (2^+2-4)7(2-^ -hi) « 4 (for 
L 3> 1), (iii) the mean shortest path length (l) increases 
logarithmically with Nl as (/) cx logiV^, and (iv) the 
clustering coefficient is zero. 

Boettcher et al. studied bond percolation on the HN- 
NP with open bond probability p [l^. In the HN- 
NP with L generations, they considered the renormal- 
ization of four parameters: Rl (the probability that 
three consecutive points a, 6, c of a chain arc connected), 
Sl (the probability of be being connected, but not a), 
Ul (the probability of ae being connected, but not b), 
and N'j^ (the probability that there are no connections 
among a, b, and c). From the renormalization group flow 
for Rl they determined the two critical probabilities as 
w 0.319445181 and p^cz ^ 0.381966. 



III. ANALYTICAL CALCULATION FOR THE 
SKELETON OF THE HN-NP 

The fractal exponent ?/'max(p) is useful to determine 
phase behavior [T^. It is defined to be ■0max(p) = 
limAT^^oo logTVf, Smax(-/VL;p). A non-percolating phase, 
a critical phase, and a percolating phase are character- 



ized by IpmaAp) = 0, < V'max(p) < 1, and IpmaAp) = 1, 

respectively. Unfortunately, it seems difficult to evaluate 
SmaxC-^Llp) directly for the HN-NP. Instead, we focus on 
a subgraph of the HN-NP and evaluate its fractal expo- 
nent. 

We extract a subgraph from the HN-NP with L gen- 
erations by removing the backbones. Because the re- 
sulting subgraph has no cycles and the number of short- 
cuts is 2^ — 2 = Nl — 3, this subgraph is composed 
of three disconnected trees. Indeed, nodes i = 0,2^~^, 
and 2^ belong to the three different trees. We call these 
the root nodes. Here the graphs isomorphic to these 
trees and having root nodes i = 0,2-^"^, and 2^ will 
be called the skeletons Ta{L), Tb{L), and Tc{L), respec- 
tively. Clearly, Ta{L) and Tc{L) are also isomorphic to 
each other, TaiL) ~ Tc{L). At L = 2, Ta{2) is composed 
of nodes and 3 and the edge between them, Th(2) is 
one isolated node i = 2, and Tc(2) is composed of nodes 
1 and 4 and the edge between them. The skeletons Ta{L), 
Tb{L), and Tc{L) for arbitrary L are given recursively as 
follows. First, we consider two subgraphs of the sets of 
nodes {0,1,- •• , 2^"!} and {2^-i, 2^-i-hl, • • • ,2-^} after 
removing the backbones from the HN-NP with L genera- 
tions. By symmetry, both subgraphs consist of the skele- 
tons T^li - 1), Tfc(L - 1), and rc(^ - 1) (FiglD^a)). Here 
the root nodes of the latter subgraph are i = 2^~^ (for 
Te(L-l)), 3x2^-2 (forTfc(L-l)),and2^ (for T,(L-1)). 
Note that the skeletons Ta{L), Tb{L), and Tc{L) are given 
by adding the two long-range edges {0,3 x 2^~2| and 
{2-^~2,2-^}, and taking into account the connection of 
node 2^-1 (FigHJb)), we have 

TaiL) ~ Te(L)~i?i(r,(L-l),rb(L-l)), (1) 
niL) ~ i?2(r,(L-l),Te(L-l)), (2) 

where the operation Ri{x, y) adds the edge between root 
nodes of the skeletons x and y, and R2{x, y) merges two 




T„(L-J) Ti(L-l) T„(L-1) T„(L-1) nil - \) T^{L - \) 
(b) 




T„(L-1) Tt(L-l) T,(L-l) T,(L-l) n{L - I) T„(L-1) 




Ta{L) n{L) T,{L) 



FIG. 2: Construction of the skeletons Ta{L), Tb{L), and 
Tc{L). Open circles represent the root nodes of the skeletons. 
See the main text for explanation. 
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root nodes of the skeletons x and y into one. 

We now calculate the mean size of a cluster including 
the root node (the root cluster size) for each skeleton. 
We denote the root cluster sizes of Tq(L)(~ Tc{L)) and 
Tb{L) by Sa{L) and si,{L), respectively. Because of the 
recursive structure ([l][2]) of the skeletons, the root cluster 
sizes Sa{L) and si,{L) also satisfy recursive relations: 

Sa{L + l) = Sa{L)+pSb{L), (3) 

Sb{L + l) = 2sa{L)-l, (4) 

where the initial conditions are Sa(2) = 1+p and S6(2) = 
1. Then, we find that 

'"^^^ " 2+ 2^+Vl + 8p 

(l + 7r+8^)^-(l-VT+8^)^ 
"'(^^-' + 2WfTW ■ 

For L > 1, wc obtain SaiL) cx iVL''''°°''"°"^^\ where 



ponent V'Sax ^^i^Llp) at a finite generation L is evalu- 
ated as 



HN-NP 



{Nl;p) 



loggmax(^L + i;p) " log Smax(^L-l ; p) 
log - log A^L-i 



(8) 

We also plot the fractal exponent ■0root°'°" (-P) of the skele- 
ton (Eq.(l7]), shown as the thick-dashed line) and p^i""^ 
and Pc2^^ (shown as vertical lines) in FigI31 

From Figiljb), we see that V'roor°"b) is actually the 



lower bound of V' 
particular, ip. 



HN-NP 



HN-NP 



(p), implying that Pd =0. In 



{Nl ; p) coincides with ip: 



skeleton 



(p) for 



p < 0.26 (except near p = 0, where finite size effects are 
not negligible). For p > 0.26, {Nl;p) is consid- 

erably greater than '(/'root"*"" ip) i and reaches unity at p = 
p^2^^. At a glance, in the large size limit, -0™^^^(A^i;p) 
seems to change continuously with p < pf^'^. However, 
we speculate that in the large size limit ^™~^^{NL-p) 
coincides with '4'ioot'^°'^{p) i^. the entire region below 
P?f^i jumps to a higher value at p = pfi^'^, increases 



monotonically again for p such that pf[^'^ < p < p^2 
and reaches unity at p = Pc2^^. The coincidence be- 



BCZ 



^T'^'ip) = iog2 (1 + v/r+8^) - 1- (7) 

We expect V'roSt"*°"(p) = '0max'°"(p) because the roots 
are hubs. In fact, we performed Monte Carlo simulations 
for the bond percolation on the skeletons. Our numer- 
ical result of "iAmax"*"" (-P) shows a good correspondence 
with Eq.© except near p = (not shown). Accord- 
ing to Eq.([7]), '0root°*°"(p) increases continuously from 
V-root^^^lO) = to V?ool°*™(l) = 1- This means that the 
subsystem consisting only of the shortcuts is in the crit- 
ical phase for all p{^ 0, 1), like the growing random tree 
[l5| . Because the HN-NP is obtained by adding the back- 
bones to the skeletons, the clusters in the skeletons be- 
come larger. Therefore, the entire system permits a criti- 
cal phase even for infinitesimal p, i.e., the non-percolating 
phase does not exist except at p = 0. 



IV. NUMERICAL CHECK 

In the previous section, we evaluated the root clus- 
ter size of the skeleton to show that its fractal exponent 
V'root'*"" (-P) takes a non-zero value for all p > 0. Because 
the skeleton is just a subgraph of the HN-NP, tpt^ot^"" {p) 
is a lower bound for the fractal exponent of the largest 
cluster of the HN-NP V'Sax"^^(p), i-e-, V'r'ool"*™^) ~ 
V'mfx*°"(p) < V'Sax"^^(p)- For bond percolation on the 
HN-NP, V'max"'^^(p) > when p > 0, implying that 
Pel = 0. To check our prediction, we performed Monte 
Carlo simulations of bond percolation on the HN-NP. 
The number of generations is L = 13, 14, • • • , 20, and 
the number of percolation trials is 100000 for each p. 

Figures (ni^a) and (b) show the results for the order pa- 
rameter m{NL]p) and the fractal exponent of the largest 
cluster {N L] p) , respectively. Here the fractal cx- 



(a) 

m 



1.0 




FIG. 3: (a) Order parameter m{NL',p) = Sui!i^{NL\p) /Nl 
and (b) fractal exponent tp^N^w)- The numbers of gener- 
ations L are 19 (black-diamond), 18 (red-circle), 17 (green- 
square), 16 (blue-triangle), 15 (open-diamond), and 14 (open- 
circle). The vertical solid line and dashed line indicate p^^^ 
and p^2^ respectively. In (b), the fractal exponent ^roSt''*°"(p) 
of the skeleton given by ((7| is shown by the thick-dashed line. 
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twccn Tpn.'^^^^ip) and V-root"'™ (p) for p < p^^'^ means 
that the partial ordering (in the sense that the largest 
cluster is 0{N") with a < 1) in this region is essentially 
governed by the shortcuts. Because Boettcher et al. [1^ 
considered renormalization of the connecting probability 
of consecutive points of the backbones, we would expect 
their first critical probability p^i^^ to be the probability 
above which the backbones become relevant. Note that 
such a jump in the fractal exponent implies a qualita- 
tive change in the criticality, and this has already been 
observed in site-bond percolation on the decorated (2,2)- 
flower (20| . In addition, our numerical result shows that 
V'niax~^^(p) reaches unity smoothly at Pc2"^. This indi- 
cates that the phase transition to the percolating phase 
is discontinuous, similarly as in [isj . 

Finally, we discuss the cluster size distribution func- 
tion, ns{p), below p^-p^- Figure HJa) shows ns{p) for 



several values of p with < p < p^^^ ■ 
phase, we expect a power law for ns(p): 



In the critical 



where 



-(p) = l + 7/Wx(p)-\ (10) 



and a corresponding scaling form: 



ns{NuP) = N-^-^^^^^^^fisN-^--^''^), (11) 



where the scaling function /(•) behaves as 



. , , I rapidly decaying func. for x 3> 1, , 
fix) ~ < _^(p) , ^ 12) 

X for X 1. 



-t(p) 



(9) 




We tested this scaling for < p < 0.26 and p^i < 
p < _Pc2"^ and obtained excellent collapses (FiglUJb)). 
We would also expect that Ug to be fat-tailed for 0.26 < 
p < p^i'^ because Ug [p) for the skeletons perfectly obeys 
Eqs.® and dTU]) via Eq.dTT]) for < p < 1 (not shown), 
and Ug is broader when we add the backbones to the 
skeletons, i.e., for the original HN-NP. 



SUMMARY 



In this paper, we have studied bond percolation on the 
HN-NP. Our results give the two critical probabilities as 
Pel = 0(< p^i^) and Pc2 = Pc2"^, implying that the sys- 
tem has only a critical phase and a percolating phase, 
and does not have a non-percolating phase for p > Q. 
As far as we know, all complex network models with a 
critical phase have only the critical phase and t he p er- 
colating phase ([sl-fisj for percolation and [ill [2ll - [26l | for 
spin systems). It will be challenging to clarify the origin 
of such universal behavior. 



FIG. 4; (a) Cluster size distribution ns{NL;p) a.t p — 0.16, 
0.26(< p^f^), and 0.36(> Pcf^), from left to right, (b) 
Scaling result for ns{NL\p) at p = 0.20(< p^x"^), 0.32, and 
0.35(> p^i^), from right to left. The numbers of generations 
L are 19 (black-diamond), 18 (red-circle), 17 (green-square), 
16 (blue-triangle), 15 (open-diamond), and 14 (open-circle). 
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